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In a quantum dot with dozens of electrons, an approximation beyond Tamm-DankofF is used 
to construct the quantum states with an additional electron-hole pair, i.e. the "excitonic" states. 
The lowest states mimic the non-interacting spectrum, but with excitation gaps renormalized by 
Coulomb interactions. At higher excitation energies, the computed density of energy levels shows an 
exponential increase with energy. In the interband absorption, we found a background level in the 
quasicontinuum of states rising linearly with the excitation energy. Above this background, there 
are distinct peaks related to single resonances or to groups of many states with small interband 
dipole moments. 

PACS numbers: 73.21.La, 78.67.Hc, 71.35.Cc 
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Semiconductor quantum dots are like anisotropic 
Thomson atoms. Most of their electronic and optical 
properties are, by now, well studiecW The typical en- 
ergy scales determining these properties are the following: 
the in-plane confinement strength, fkv, which also defines 
the unit of length, the characteristic Coulomb energy, (3, 
and the characteristic energy along the symmetry axis, 
E z . The ratio (3/ftw is an indicator of the role played 
by Coulomb interactions in the dot. In the strong con- 
finement regime, (3/huj <g; 1, which we encounter in self- 
assembled structures, quantum dots are like real atoms, 
with correlation energies reaching only a few percents of 
the total energy. On the other hand, for (3 j 'hui ^> 1, corre- 
lations are strong, and there could be even premonitions 
of Wigner crystalization for large enough dots. 

The excited states of a quantum dot, as in most quan- 
tum systems, subscribe to a general picture of a few low- 
lying "oscillator" states followed by a quasi-continuum 
of excitations. Typically, the low lying states determine 
most of the optical properties. But there are situations in 
which higher excitations play a key role. Raman scatter- 
ing or PLE experiments with excitation energies 40 — 60 
meV above band gap, for instance, sense high interband 
excitations. On the other hand, the density of energy 
levels in the quasi-continuum obeys a simple exponential 
rule with excitation energy, known in Nuclear Physics as 
the Constant Temperature Approximation (CTA^P. I n 
paper [3], we corroborated this picture for the intraband 
excitations of few-electron quantum dots in a magnetic 
field. 

In the present paper, we focus on the interband exci- 
tations. A medium-size dot, charged with 42 electrons, 
under intermediate confinement regime (P/hiu w 0.6) is 
studied, a situation typical of etched dots. As it will be 
shown below, the obtained low-lying spectrum is similar 
to the spectrum of the non-interacting model, but with 
values renormalized by Coulomb interactions. The quasi- 
continuum of states follows the CTA, as expected. 

For a dot with 42 electrons, no exact diagonalization 
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FIG. 1: (Color online) The pp-TDA and ppph contributions 
to the wave function of the excitonic states. 



calculations are possible and we should resort to approxi- 
mations. We will use a kind of Quantum Chemistry Con- 
figuration Interaction Method (CIp- in which the start- 
ing point is the Hartree-Fock (HF) solution for the 42- 
electron problem. Excitations over the HF solution con- 
form the basis functions in which our Hamiltonian is to be 
diagonalized. The simplest set of excitations, one addi- 
tional electron in the conduction band (CB) and one hole 
in the valence band (VB), is usually called pp-Tamm- 
Dankoff approximation (pp-TDA) in Nuclear Physics 5 . 
The valence-band hole is treated as a quasiparticle with 
its own dispersion relation. Excitonic states are, thus, 
states with two additional "particles" above the HF so- 
lution. The next set in increasing complexity includes 
excitations with an additional electron above the Fermi 
level and a hole below the Fermi level in the CB. We will 
call it the ppph contribution. A schematic representa- 
tion of the pp-TDA and ppph contributions to the wave 
function of the excitonic states is given in Fig. [T] We 
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will truncate the basis of functions at this stage. 

The reason for including ncxt-to-lcading order excita- 
tions in the basis functions, i.e. going beyond pp-TDA, is 
the following. It is well known that the pp-TDA underes- 
timates the density of energy levels already for relatively 
low excitation energies. A correct level density could 
be crucial in the description of Raman scattering events 
with laser excitation energy well above band gap, which is 
our next goal. A particularly interesting question is why 
the so called outgoing resonances, presumably a result of 
higher order Raman processes^, give stronger contribu- 
tion to the Raman spectrum than the usual second-order 
Raman amplitude. 



II. FORMALISM 

As mentioned above, the starting point in our for- 
malism is the HF solution for the 42-electron problem. 
The HF equations and the way we solve them are de- 
scribed elsewhere^'. We use the following set of pa- 
rameters: m e = 0.067 mo, (a GaAs quantum dot), 
hw e — 12 meV (the confinement energy for electrons), 
and /3 = 7.07 meV. The latter value comes from the ex- 
pression: j3 — 0.6 e 2 /(47reo«:Z), where the relative dielec- 
tric constant is k = 12.8, I is the oscillator length, and 
the factor 0.6 takes account of the quasibidimcnsionality 
of the structure. Only one electron sub-band in the z- 

(e) 

direction is considered, and the value of E\ is used to 
control the effective band gap of the nanostructure. 

The HF states for the hole in the VB are found from the 
Kohn-Luttingcr Hamiltonian, in which a term accounting 
for the background potential created by the 42 electrons 
is included^. 

For the intcrband excitations of the dot, i.e. many- 
particle wave functions with 43 electrons and one hole, 
we write the ansatz, 



F = l a + l T — rrij 

= l a + l p - l\ + l T - rrij , (2) 



— Su + Sp S\. (3) 

Here I are (orbital) angular momentum, rrij - (band) 
hole angular momentum, and S - (electron) spin quan- 
tum numbers. Of course, these magnitudes correspond 
to projections along the z-axis. 

The coefficients V and Z (the wave functions) and the 
excitation energies with respect to the HF state are com- 
puted from the following eigenvalue problem, in close 
analogy with the TDAE1 

[H,Q^]\HF) = AE int Q^\HF), (4) 
where H is the many-particle electron-hole Hamiltonian 

H = T e + T h + V ee + V eh . (5) 

The T are one-body operators, and V refer to Coulomb 
interactions. Let us stress that electron-hole exchange, 
usually of the order of peV in these structures, is ne- 
glected. The explicit matrix form of Eq. Q is presented 
in Appendix A for completeness. We use energy cutoffs of 
90 meV to reduce the dimensions of the resulting matri- 
ces to around 12000. Eigenvalues and eigenvectors (the 
wave functions) are obtained by numerical diagonaliza- 
tion. 



III. THE LOWEST EXCITONIC STATES 
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a and p are electronic states above the Fermi level, and A 
is a state below the Fermi level of electrons. We restrict 
the second sum in Eq. (JT|) to a > p in order to avoid 
overcounting. This second sum is the term beyond pp- 
TDA, not included in previous computations of interband 
excitations^'. 

Due to the peculiarities of Coulomb interactions, the 
Kohn-Luttinger Hamiltonian, and the central confine- 
ment potential, the interband excitations are character- 
ized by two quantum numbers, F and S z . It means that 
the single-particle states entering the first and second 
sums of Eq. should satisfy, respectively, the first and 
second rows of the following equalities: 



The solution of Eq. |4 ) provides us the wave functions 
and energies of the interband excitations of the dot (the 
excitonic states). In this section we focus on the lowest 
states. 

A sample of results is shown in Fig. [2ja). Energies 
are measured with respect to Ehf + Ez - States with 
quantum numbers F = —3/2, S z = 1/2 are represented 
in the figure. With increasing energy we observe first 
a pair of states followed by a group of four states, and 
then a quasi-continuum of excitations. This sequence can 
be understood in terms of the non-interacting picture of 
electrons and holes. 

Our model dot with 42 electrons has 6 filled 2D os- 
cillator shells. We should add an electron-hole pair to 
construct an interband excitation. The electron should 
then be added to the 7th shell. The possible angular 
momentum values are l e — —6, —4, . . . , 6. The hole or- 
bital state, on his hand, can be the first oscillator state 
with lh = 0. As we are talking about F = —3/2 states, 
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FIG. 2: (a) The low-lying interband excitations in the model 
quantum dot. States with F = —3/2, S z = 1/2 are drawn. 
Results from the non-interacting picture are also shown for 
comparison, (b) Increase in the density of energy levels due 
to the ppph component. 
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FIG. 3: (Color online) (a) Increase in the number of levels 
as predicted by the CTA for states with 9 meV < AE < 15 
meV. (b) Probability distribution of level spacing for excitonic 
states in the interval 9 meV < AE < 12 meV. 



only two combinations remain: l e — 0, mj — 3/2 and 
l e = —2, mj = —1/2. That is, one heavy and one ligth 
hole exciton that, because of the assumed 25 nm width 
of the dot in the z-direction, are very close in energy (1.5 
meV). The next set of 4 states are related to hole exci- 
tations, i.e. the hole occupies one oscillator state of the 
2nd shell. We choose the confinement frequency for holes 
in such a way that the oscillator length is unique, mean- 
ing that electrons and holes are confined in the same 
spatial region. For heavy holes, for example, it means 
that tOhh = rn e ui e /mhh- As a result, excitations energies 
for holes are smaller, of around 5 meV. Next, we reach 
a group, more numerous, of states built up from one-huj 
electron excitations (12 meV) or two-hw hole excitations. 

Following this non-interacting picture, we can assign 
the first two levels to heavy-hole and ligth-hole excitons, 
the next four to hole excitations to the 2nd oscillator 
shell, and we can relate the quasi-continuum of states to 
electron excitations to the 8th shell and hole excitations 
to the third shell. The effects of Coulomb interactions are 
apparent in Fig. [2|a): energy values are pushed down, 



energy gaps are reduced and a quasi-continuum of levels 
appear at excitation energies even below Htu e = 12 meV. 

In Fig. |2jt>) we illustrate the main effects of the ppph 
component on the energy levels. As compared with the 
pp-TDA, we get small variations of the ground state en- 
ergy (of around 3 meV) and the gap to the first quartet 
of states. But, for excitation energies above 10 meV a 
significant increase in the level density is observed. In 
the next section a more detailed study of the density of 
energy levels is carried on. 



IV. THE CONSTANT TEMPERATURE 
APPROXIMATION (CTA) 

It can be verified that the CTA is a good approxima- 
tion for the density of energy levels at relatively small 
excitation energies, starting from the quasi-continuum of 
states. 

We show in Fig. [3ja) the excitonic sector with quan- 
tum numbers F = —3/2, S z = 1/2. Excitation energies, 
measured with respect to the lowest state in this sector, 
run in the interval 0-15 meV. A perfect exponential be- 



haviour is apparent in the regions 9 meV < AE < 12 
meV and 12 meV < AE < 15 meV. In these intervals, 
the number of states can be fitted to: 



N = N exp 



(¥)• 



(6) 



where: N = 4.472 x 1(T 3 , 6 = 1.268 meV for 9 meV 

< AE < 12 meV, and N = 1.145, 6 = 3.183 meV for 
12 meV < AE < 15 meV. 

We guess that a simple qualitative relation should exist 
between the temperature parameter, (the slope), and 
the system parameters. Indeed, in a scaled formulation 9 
should be proportional to Hu> e , simply because of di- 
mensional arguments. The number of electrons, N e , and 
the scaled Coulomb strength, ft = ft/(Huj e ) should appear 
in the combination iVe 4 /3P We expect the following de- 
pendence: ~ huj e /(Ne^ 4 ft) 2 ^ 3 . A verification of the 
above dependence requires extensive calculations and is 
delayed for a future work. 

The slope discontinuity at AE sa 12 meV is a remi- 
niscence of shell structure. Indeed, if there were a gap 
between perfect shells, then the curve would be flat im- 
mediatly after AE sw 12 meV. Coulomb interaction is 
responsible for filling the "gap" region, 12 meV < AE < 
15 meV, with levels. 

The interpretation of levels in the interval 9 meV 

< AE < 12 meV as a "distorted shell" is also reinforced 
by the distribution of spacing between neighbourgh lev- 
els, as shown in Fig. [3jb). The probability density vs 
spacing, normalized to the mean level spacing, which in 
the present situation is 0.067 meV, is shown to follow 
a Poissonian dependence, indicating a regular region of 
phase spaced This statement proves that the CTA is 
not necessarily related to chaos in the excitation spec- 
trum. 

In conclusion, the quasi-continuum of excitonic states 
for 9 meV < AE < 12 meV w Hu> e exhibits regular or 
quasi-integrable behaviour and is well described by the 
CTA. Above this interval, an abrupt change in the tem- 
perature parameter is reminiscent of shell strucure in the 
model. 



V. INTERBAND ABSORPTION 

Interband absorption at normal incidence is charac- 
terized by the creation of electron-hole pairs mainly in 
monopole states, l e + lh = 0. This means low values of 
\F\. " 

The transition matrix element for absorption is given 

by 
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FIG. 4: (Color online) (a) Interband absorption in our model 
dot and the density of energy levels, (b) Contribution to 
the absorption by the set of states with quantum numbers 
F = -3/2, Sz = 1/2. 



where the initial state is the HF state for the 42-electron 
system, and is the light-matter coupling Hamilto- 

nian corresponding to the absorption of a photon and cre- 
ation of a pair. An explicit expresion for is given in 
Refs. [8 10J. The matrix elements b aT are explicitly given 
in Ref. [7], where they are termed band-orbital factors 
(see Eq. (30) of that paper). Notice that only the first 
part of the wave function ([l]) enters the matrix element 
([7]) . This is the result of treating the ground state in the 
Hartree-Fock approximation, lplh corrections would not 
modify Eq. (|7| because the Hartree-Fock state is orthog- 
onal to lplh corrections^. 

We show in Fig. |4][a) the computed zero-temperature 
interband absorption in our model quantum dot. We 
consider normal incidence and linear polarization. The 
absorption intensity has contributions from four sets of 
states, that is (F, S z )=(-3/2,l/2), (3/2,-1/2), (-1/2,-1/2) 
and (1/2,1/2). A Lorentzian broadening of the spectral 
lines is assumed. The energy widths, Ti nt , are choosen 
in the following way: Tint — 0-5 meV for AE < 35 meV, 
and T int = 2.0 for AE > 35 meV, where AE is the 
energy difference with respect to the first excitonic state, 
and the value 35 meV is the LO phonon energy in GaAs. 
Our ansatz for Tint expresses the fact that the opening 
of the channel for the emission of a LO phonon broadens 
the excitonic levels. 
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The absorption spectrum has the usual form. First, a 
few isolated peaks, and then a quasicontinuum. We have 
superimposed in Fig. |4^a) a histogram with the den- 
sity of energy levels, in logarithmic scale. To construct 
the histogram, we used an energy width equal to Ti nt , 
that is 0.5 meV. The decrease of the density of states for 
AE > 45 meV is an artifact related to the truncation 
of Hilbert space. It becomes evident that in the qua- 
sicontinuum there is a background absorption intensity 
roughly proportional to the logarithm of the level density, 
i.e. to the excitation energy AE (because of the CTA). 
We notice that this is the dependence that follows for ab- 
sorption in the independent-particle picture for electrons 
in a harmonic potential. Above this background level, 
we observe distinct peaks. 

The contribution to the interband absorption by the 
set of states with quantum numbers F = —3/2, S z = 1/2 
is sketched in Fig. (4^b). The comparison with FigQa) 
reveals that there are peaks associated to one (or a few) 
excitonic states with strong interband dipolc moments, 
but there are also peaks due to the contribution of many 
states with relatively small dipole moments. The lat- 
ter are better understood as broad resonances, but they 
should be distinguished from the background intensity 
created by the "single-pair" excitations. 

VI. CONCLUDING REMARKS 

In the present paper, we have computed the excitonic 
states in a relatively large quantum dot by using a Con- 
figuration Interation scheme. This is an ab-initio pro- 
cedure, fully variational. The only approximations in it 
are the effective mass description of electrons, the Kohn- 
Luttinger description of holes in the valence band and 
the harmonic model for the dot. To the best of our 
knowledge, this is the first calculation of such kind in 
the literature. 

For the sake of simplicity in the presentation, we used a 
closed-shell situation, 42 electrons, in which the ground 
state is characterized by quantum numbers L to tai = 0, 
and 5 = 0. However, our scheme works properly also for 
open-shell dots and, with appropiate 2p2h corrections, 
provides quantum numbers for the ground state in accor- 
dance with Hund's rules^. Comparison with exact calcu- 
lations for few-electron dots shows excelent agreement!^. 



The main results of the paper are the following: (a) 
A picture for the low-lying states in which the role of 
Coulomb interactions is made evident, (b) A universal 
parametrization of the density of energy levels for exci- 
tation energies greater than 1 hu, and (c) A picture for 
the interband absorption made up from peaks of differ- 
ent nature and a background level in the quasicontinuum 
rising linearly with energy. 

The work may be continued along different lines. A 
natural extension is the description of biexcitonic states. 
The verification of universality of the CTA in the level 
density of quantum systems is also worth trying. Finally, 
the obtained excitonic states can be used in the computa- 
tion of Raman cross sections for laser excitation energies 
well above band gap. Reasearch along these directions is 
currently in progress. 
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APPENDIX A: EQ. Q IN MATRIX FORM 

Eq. ([l]) for the operator makes evident the basis 
functions used in the construction of the excitonic states. 
They are of two kinds: 

etht\HF), (Al) 

and 

eU$e x ht\HF). (A2) 

If Eq. Q is projected onto these basis functions, we 
get the following system of equations: 



J2{HF\[h T ,e„,, [H,^h\]]\HF)V aT + £ (HF\[h T ,e a >, [H,eUlexht]}\HF)Z rTpXT 



apXr 



AE mt J2 {HF\[h T , e„, , 4 hi } | HF) V ar , (A3) 



Y J (HF\[h T ,e\,e p ,e a ,,[H,elh\]]\HF)V (Tr + £ {RF\\h T ,e\,e p ,e a , ,\H,e\e\e x h\\\\UF)Z apXr 

o"r apXr 
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AE mt Y / (HF\[h T ,e[,e p ,e a> ,eUlexhl}\HF)Z apXT , (A4) 

crpXr 



where AE int is the excitation energy, measured with re- The matrix elements entering Eq. (A5l are computed by 
spect to the ground state of the iV e -clcctron system, and means of the following expressions: 
H is the quantum dot Hamiltonian. These equations can 
be rewritten in a more compact form: 

& b)( V z)= ae ^( V zI ^ 



(ar\A\a'T') = (e ff + e T )5 aa ,6 TT , - /3(aY|~K) , (A6) 

Teh 



{apXr\B\a' p'X'r') = (e a + e p + e T - e\)5 aa ,5 pp >5 TTl 5x\' 



13 ({p'a'\-\pa)8 TTl 5 xx , + {Xp'\-\pX')5 TT ,5 aal + (X*'\-\aX')6 rT ,5 pp , 
V r e r e r e 

(Xp'\-\X^)6 TT ,8 pa , + (Xa'\±\X r p~)8 TT/ 5 ap ,) - p ( {u't'\^-\ot) 5 pp ,6 xx , 
(p'r'\ — |ot> Vfcy ~ (At 7 | — |AV)V«^ - {o , T J \ — \pr)5 trf ,6x X , 

i eh f*eh 'eh 

(p't'\ — |pt)<W<5aA (A7) 

Teh J 



(apXT\C\a'T') = p ( (ff'AI-lpff)^ - (Ar'| — |or) 5 pa , + (Ar'| — |pr)<W ) , ( AM 

• e '"eh, ^eh, 



(ar|C>VAV) = /? ( ( P V|-|aA')^ - (aV| — |A't)V + <A'| — |A'r)<W ) , (A9) 

' e '"eh, ^*e/i 



i 

where we used the notation (ij\\kl) — (ij\\kl) — of the Hartree-Fock state i, and (3 is the Coulomb inter- 

for the Coulomb matrix elements, ej denotes the energy action strength. 
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